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Abstract: We develop a framework of calculating entanglement entropy for non-conformal 
field theories with the use of the dilaton effective action. To illustrate it, we locate a theory 
on a cylinder M x § 2 and compute entanglement entropy of a cap-like region perturbatively 
with respect to the mass for a free massive scalar held. A renormalized entanglement entropy 
(REE) is proposed to regularize the ultraviolet divergence on the cylinder. We find that the 
REE decreases monotonically both in the small and large mass regions as the mass increases. 
We confirm all of these behaviors by the numerical calculations, which further shows the 
monotonic decrease of the REE in the entire renormalization group flow. 
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1 Introduction 

Quantum aspects of field theories manifest themselves in various forms including anoma¬ 
lies, renormalization group (RG) flow, and strong dynamics on the long distance. Some of 
them can be uncovered through correlation functions of local operators that are only sen¬ 
sitive to the short-scale structure, while non-local observables like Wilson loops are needed 
to detect long-range correlations in quantum field theories (QFTs). Among such non-local 
quantities is entanglement entropy (EE) that diagnoses quantum entanglement relevant at 
low energy where thermal (classical) fluctuations are suppressed enough. In critical phe¬ 
nomena, EE shows universal and characteristic features that reflect the central charges of 
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conformal field theories [1-5] and the degeneracy of the ground state of topological QFTs 
[6-8]. It also plays an important role as a measure of degrees of freedom under RG flows and 
confinement/deconfinement transitions [4, 5, 9-16]. 

Let us assume that the total system is in its ground state | fl). We denote the subsystem 
for which we want to compute the EE by A. The entangling surface is the boundary of A 
which we denote by £ = dA = dA. The density matrix describing the total system is given 
by 

p vac = \n)(n\. (i.i) 

If the Hilbert space of the theory is assumed to factorize as T~i = T~La < 8 > T~La, one can define 
the reduced density matrix of A by 


PA = tl'H 4 Pvac ■ (1-2) 

The EE of the subsystem A is then defined as the Von Neumann entropy associated with the 
reduced density matrix pa, 

S A = — tr Ha Pa^pa • (1.3) 

In QFTs, EE is usually calculated by the replica trick. Applying the replica trick, EE 
can be computed using the formula [1, 3] 

S A = lirn d n (F n - nFi) . (1.4) 

n— S>1 

Here F n = — In Z[A4 n ] is the free energy of a theory on the n-fold cover Ai n of the Euclidean 
space-time on which the theory lives. Let us describe this in some detail. 

The entangling surface E is a codinrension-two hypersurface in the (Euclidean) space- 
time. Parameterizing the transverse two-dimensional space by the polar coordinates (p, r) so 
that £ = {x\p = 0}, the n-fold cover M. n is obtained by extending the periodicity of the 
angular coordinate r from 2n to 27rn. For n/ 1, the n-fold cover A4 n has a conical singularity 
at £ of the angular excess 27r(n — 1). This singularity prevents us from calculating the exact 
free energy F n of QFTs except for free field theories [17, 18]. In practice it is not easy to work 
on such singular spaces in particular when the translational symmetry along t is broken in 
the transverse space. A common practice is to regulate the conical singularity by replacing 
the conifold with a smooth manifold with a regulator parameter and remove the regulator 
at the end of the calculation. This procedure is ambiguous because there are more than one 
methods of regularization and it is not clear whether all of them give the same answer for 
EE. Instead one could work on the singular conifold without regularization. 

It is, however, not so obvious how to implement conformal perturbation theory on a cone. 
Proper determination of boundary conditions at the singularity is difficult. In this paper we 
will attack this problem in a different way. We will call it Conformal Decompactification. 
The idea is to perform a conformal transformation of the replica space such that the resulting 
space has no conical singularity and calculate the free energy on the non-singular geometry. 
The only constraint on the conformal factor is that it should blow up along the entangling 
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surface where the conical singularity is located. For non-conformal field theories the two free 
energies are made equal by switching on a background dilaton field on the non-singular space. 

In section 2, we will describe the conformal decompactification in a relevant perturbation 
of conformal field theories (CFTs) and provide a framework for the perturbative calculation 
of the EE. In general, the perturbation theory of EE on non-compact spaces is plagued by 
infrared (IR) divergences. In order to avoid such IR divergences, we will locate the theory on 
a compact space such as a cylinder M x with the metric 

ds 2 = — dt 2 + R 2 (d6 2 + sin 2 9 dO^_ 2 ) , (1.5) 

and consider the EE of a cap-like region A = {x|t = O,O<$<0o}- This type of EE has 
been extensively studied in recent literatures [19-22], but it remains unclear how to calculate 
the entropy and how it behaves away from conformal fixed points. We wish to explore aspects 
of the EE on the cylinder (1.5) in non-CFTs. 

For brevity, we will confine ourselves to the case in d = 3 dimensions, where the entropy 
on the cylinder (1.5) has the area law divergence proportional to RsuiOq that needs to be 
regularized to extract a scheme-independent part from EE. The situation is the same for EE 
of a disk of radius i?disk on a flat space and we briefly review the known facts so as to extend 
them to the case we are dealing with. 

A simple and pragmatic regularization of the UV divergence of EE for a disk is proposed 
by Liu and Mezei [10] as the renormalized entanglement entropy (REE): 

-T(Rdisk) = (Rdisk3i? disk - l)S'A(Rdisk) ■ (1-6) 

Among the other various regularizations, the REE turns out to have particularly important 
properties as a measure of degrees of freedom in QFT: 

1. It is monotonic T' < 0 and the equality holds only for a theory being conformal [11]. 

2. The value of T for CFT 3 is the same as the free energy on a three-sphere, E = — log 3 

[23]- 

Upon the identification of the size of the disk with the inverse of the renormalization scale, the 
REE monotonically decreases under any RG flow, proving the E-theorem that states Euv > 
Fir for the free energies at the UV and IR conformal fixed points [15, 16, 24, 25]. We will 
call a function satisfying these properties an E-function. An E-function is reminiscent of the 
Zamolodchikov’s c-function [26] that also satisfies stationarity in addition to the counterparts 
of the two properties . 1 In our case, the stationarity of an E-function can be stated in the 
following way: 

3. The derivative of T with respect to a coupling constant g of a relevant operator vanishes 
at the conformal fixed points of RG flows, <9 9 E|cft = 0. 

lr The second property of the Zamolodchikov’s c-function is that the value at the fixed points is the same as 
the central charge of the CFT 2 . 
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The third property, not necessary for the proof of the E-theorem, can be favorable for finding 
the fixed points of RG via an E-function. The stationarity of the REE (1.6), however, is 
doubted for a massive free scalar theory based on the numerical calculation [27, 28] where it 
is speculated that the non-stationarity of the REE originates from the IR divergence of EE. 

One of the aims of this paper is to seek an E-function in the Zamolodchikov’s sense, i.e. 
, satisfying all the aforementioned three properties, through the EE on a cylinder which is 
free from the IR divergence. 

There are two natural generalizations of the REE on the cylinder as will be introduced 
in section 3. One is the Liu-Mezei type, called Elm, which takes the same form as (1.6) with 
Rdisk replaced with the size of the sphere R while fixing the angle 6 q. On the other hand, 
one can define another function, called E<c, by replacing i?disk with sin$o hi (1.6) while fixing 
R, to remove the UV divergence. We will show both the REEs on the cylinder are finite, 
coincide with the three-sphere free energy F at RG fixed points, and reduce to the original 
REE (1.6) in the flat space limit R —> oo and do > 0 while keeping R6q = Rdisk held fixed. 
The REEs are weak E-functions in a sense that their values at the UV fixed point are larger 
or equal to those at the IR fixed point as a consequence of the E-theorem. 

Reminding the situation for the REE on the flat space, a free massive scalar field is a 
good testing ground for examining the stationarity of the REEs on the cylinder. As the 
entropy is dimensionless, Sa is a function of dimensionless parameters do and mR for a scalar 
field of mass m. In section 4, we carry out the perturbative expansions of SA{do,mR) in 
the small and large mR limits respectively. In the small mR limit, we apply the conformal 
decompactification developed in section 2 to the cylinder and retrieve the leading contribution 
of order (mR) 2 to the entropy. In the other limit, we slightly generalize the method of [29] 
that relates the order 1 / (mR) term in the large mR expansion with the logarithmic divergence 
in the EE of a four-dimensional free massless scalar field. In both regimes, we will show both 
of our REEs, Elm and Ec, decrease monotonically as mR becomes large for fixed do and as 
#o increases from 0 to 7 t/2 for fixed mR. Meanwhile, we find Ec is stationary at the UV fixed 
point (mR = 0) , being a good candidate for an E- function, while Elm still is not stationary 
even on the cylinder. 

The limiting behaviors of Ec are confirmed by the numerical calculation which further 
demonstrates the monotonicity under the entire RG flow. The numerical calculation is based 
on the real time approach [30, 31] which does not rely on the replica trick. Thus it supports 
the validities of the conformal decompactification and the large mass expansion we use in our 
analytic calculations. 

Put it all together, the examination of the EE of a free massive scalar field provides us a 
non-trivial evidence for Ec being a strong E-function that monotonically decreases under any 
RG flow. It satisfies the three properties at least for a free massive scalar field, and could be 
an E-function in the Zamolodchikov’s sense. It would be intriguing to extend the proof of the 
E-theorem [11, 32] to the cylinder, which might prove the monotonicity of our REEs. On the 
other hand, we were not able to verify the monotonicity of Elm even numerically. Moreover, 
we are led to a surprising discrepancy between the numerical result and the analytic large 
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mass expansion that arises only for J- lm- We will give the details of the numerical algorithm 
and the conjectured form of the large mass expansion fitted from the result in section 5 where 
possible sources of the inconsistency are also discussed. 

2 Conformal decompactification 

We illustrate the conformal decompactification by taking our space-time to be (2 + 1)- 
dimensional Minkowski space and the subsystem to be a disk of radius R on some spatial 
plane. The entangling surface is a circle of radius R. Since we are interested in a static 
situation we can as well do everything in the Euclidean space M 3 . In a coordinate system 
adapted to the disk the metric of M 3 can be written as 

ds ^ 3 = dp 2 + p 2 dr 2 + (R + p cos T) 2 d(j) 2 . (2.1) 

The coordinates p and r are the radial and angular coordinates, respectively, on the plane 
transverse to the circular entangling surface located at p = 0. The angular coordinate (j) along 
the entangling circle has periodicity 27 t. 

The replica trick amounts to changing the periodicity of r from 27r to 27rn, where n is a 
positive integer. If n ^ 1, then there is a conical singularity at p = 0 with the angular excess 
27r(n — 1). In this example, we have a conical singularity because the transverse r-circle is 
shrinking to zero size along the entangling surface. One way of getting a regular space is 
to perform a conformal transformation which blows up along the entangling surface. In the 
resulting conformally transformed space the transverse circle will be non-contractible with 
coordinate periodicity 2irn and there will be no conical singularity. We can read off the free 
energy on the conifold from the free energy on the conformally related smooth manifold if we 
know how the free energies are related . 2 

Let us see in the particular case of the disk how this conformal transformation can be 
done. Following [23, 34] we write the metric of M 3 in the cylindrical polar coordinates as 

ds^ 3 = dt 2 E + dr 2 + r 2 d(j) 2 , (2.2) 

where t E is the Euclidean time and the entangling circle is located at ( t E ,r ) = (0,-R). This 
is a different coordinate system from that was used in (2.1). If we now define 

u = r + it E , E = u + ir E , (2-3) 

and make the following coordinate transformation 

_v R — CU 

e - 1 =- 

R + u 

the metric of M 3 can be written as 

R 2 (dr E + du 2 + sinh 2 u dcf) 2 ) = e 2a ds \ 3 , (2.5) 

2 The reader should note that this method is similar in spirit to the method of [33]. 


(2.4) 
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where the conformal factor is given by 


e 


<J 


2 R 2 

| R 2 -uj 2 


( 2 . 6 ) 


The metric on the left hand side is the metric of two-dimensional hyperbolic space H 2 
times a circle parameterized by te with periodicity 27 t. If we extend the periodicity of te from 
27T to 27rn, then the metric of H 2 x S 1 becomes conformal to the replica geometry where there 
is a conical singularity along the entangling circle located at uj = R of angular excess 27 r ( n — 1 ). 
The conformal factor is independent of n and is a periodic function of te with periodicity 27r. 
Hence the replica geometry corresponding to a disk can be conformally mapped to H 2 x S 1 
with S 1 periodicity 27 xn and the conformal factor blows up along the entangling circle at 
uj = R. The S 1 factor is non-contractible in the resulting geometry and this S 1 is precisely 
the image of the contractible circle in the plane transverse to the entangling surface located 
at uj = R, which is used to perform the replica trick. 

Although we have explained the procedure for a disk, this method is general because 
every entangling surface is locally the same. There is a contractible circle in the plane trans¬ 
verse to the entangling surface and we can make it non-contractible by doing a conformal 
transformation which blows up along the entangling surface. The resulting geometry will 
have a non-contractible circle and will be smooth. For example, if we take as our subsystem a 
two-dimensional region bounded by the curve, r = /(</>), where /(</>) is a periodic single-valued 
function of q i, then we can choose the conformal factor to be of the form: 


' |/ 2 (</>)- u , 2 


(2.7) 


The coordinates 4> and uj are as defined before. The conformal factor diverges along the curve 
r = f(4>) located on the tE = 0 plane. The resulting geometry obtained by multiplying the 
M 3 metric by this conformal factor is no longer H 2 x S 1 but it is smooth and has a non- 
contractible circle. This circle has periodicity 27T and every field theory quantity is periodic 
under 2ir shift along the circle. The replica trick in this geometry amounts to making the 
periodicity of the circle 2nn. This does not produce any singularity because n is a positive 
integer and things are periodic along the circle with periodicity 27r. 


2.1 Relating the free energies on conformally related spaces 

For a conformal field theory the free energies on the replica space and the conformally trans¬ 
formed space are the same modulo the conformal anomaly if the space-time dimension is even. 
The conformal anomaly part can be computed by standard methods. 

For a general non-conformal field theory the free energies on the two spaces are not equal, 
but they can be related if we introduce a background dilaton field, which we denote by r(x). 
Let us denote the metric of the replica space by g^^x). 

If Z[g fiu (x), r(x)] is the partition function of the Euclidean theory in the presence of the 
background metric g^,i,(x) and dilaton held r(x), then it satisfies the following transformation 
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rule [35-38] 


Z[e 2a(x) g^(x),T(x) + a(x)} = C Z[g^(x),T(x)] , (2.8) 

where C is completely determined by the conformal anomaly of the ultraviolet (UV) CFT 
and does not depend on the mass parameters of the theory. In particular C = 1 in odd 
dimensions due to the absence of conformal anomaly. The free energy defined as F = — In Z 
satisfies the relation: 

F[e 2rr(x) g r(x) + <r(x)] = F\g^, r(x)] , (2.9) 

where we have neglected In C, because we are only interested in the part of the EE generated 
by the massive deformation. In odd dimensions this factor is identically zero and this equality 
is exact. In even dimensions this anomaly part gives rise to local terms in the dilaton effective 
action some of which are uniquely determined by the trace anomaly matching. These local 
terms in the dilaton effective action give the logarithmically divergent universal terms in the 
entanglement entropy which were computed by using this technique in [39, 40]. 

Now the equality (2.9) holds for any functional form of the dilaton field t(x) and we can 
also write 

F[e 2 °Wg^, <r(x)] = F[g^ r(x) = 0] . (2.10) 

The right hand side represents the free energy on the replica space in the absence of the dilaton 
field, which is precisely what we want to compute, and the left hand side represents the free 
energy on the conformally related non-singular space but in the presence of a background 
dilaton field which is equal to the conformal factor <r(x). We will use this relation to compute 
the EE by conformally mapping the problem to a non-singular space. 

2.2 Deformed CFT coupled to dilaton 

Let us consider a UV CFT in d dimensions deformed by some (marginally) relevant operator 
Oa of dimension A. The action can be written as 

I = Luvcft + g d ~ A J d d x yfg A(/r) Oa , ( 2 - 11 ) 

where A(/x) is the dimensionless renormalized coupling constant at renormalization scale /r. 
It is determined by the beta function equation 

^A = /3(A). ( 2 . 12 ) 

After conformally transforming to the non-singular space the free energy has to be cal¬ 
culated on the new geometry with dilaton turned on. The dilaton held couples to the held 
theory as [36] 3 

I = Ajvcft + n d ~ A j d d x \(ge T ) Oa ■ (2-13) 

3 See appendix A for an example of the action on a conformally transformed manifold. 
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If we set t = a where g^ u = e 2a g^-, we get the following action 

I = Ajvcft + /r d_A J d d x y / ^A(/re fT ) Oa , (2.14) 

which needs to be used for the calculation of free energy on the conformally transformed 
manifold. To summarize, in our prescription, one needs to calculate the dilaton (with t = a) 
effective action on the conformally related non-singular manifold to compute the EE. 

We would like to emphasize that we have not used any perturbation theory to arrive at 
this prescription. Thus it can be used even if conformal perturbation theory breaks down. 


3 Renormalized entanglement entropy on cylinder 

We will consider a theory on a cylinder M x with the metric given by (1.5) and divide 
the S d_1 by a codimension-two hypersurface E at t = 0 and 0 = Oq to a subsystem A within 
0 < 9 < 9 0 and its compliment A within 9$ < 9 < tt. The angle 9q can be restricted to be 
0 < 9q < -k/2 for the entropy is symmetric with respect to the exchange of A and A when we 
concentrate only on the vacuum state of the theory. 

Employing the replica trick, one can calculate the EE with the partition function on the 
n-fold cover of the Euclidean space of (1.5) 

ds 2 = dt 2 E + R 2 [d9 2 + sin 2 0dfl^_ 2 ) , (3.1) 


that has a surplus angle 27r(n — 1) around E. To make it transparent, we use the coordinate 
transformation 


twh.it e / R) 
tan 9 


sin #o sinr 
cosh u + cos 9q cos t 
sin 0o sinfui 
cos 0 q cosh u + cos r 


(3.2) 


with 0 < u < oo and 0 < r < 27 t for n = 1. The resulting metric becomes [23] 


e 2a ds 2 = R 2 [dr 2 + du 2 + sinh 2 u dQ^_ 2 ] , 

e - 2 a = _ sin 2 0 O _ (3-3) 

(cos r + cos 0o cosh u) 2 + sin 2 0o sinh 2 u 

The n-fold cover is given by the metric (3.3) with the period r ~ T-\-2irn, which is conformally 
equivalent to S 1 x We will denote the conformally equivalent manifold as x Hk^ 1 . 

The entangling surface E located at tE = 0 and 0 = 0o in the original coordinates (3.1) 
is mapped to u = oo in the new coordinates (3.3) where the conformal factor e 2a blows up. 
Note that the S 1 factor along r is non-contractible in the resulting geometry. This 8 1 is the 
image of the contractible circle in the plane transverse to the entangling surface at = 0 
and 0 = 0q, which is used to perform the replica trick. 





EE is always accompanied by UV divergences in QFT. The leading part is well-known 
as the area law term diverging as l/e rf ~ 2 in d dimensions for the UV cutoff e <C 1. For this 
reason, the bare entropy is scheme-dependent and needs to be renormalized so as to be free 
from the UV divergences. One possible regularization is to renormalize the divergences to 
parameters in the background gravity theory such as the Newton and cosmological constants 
as is usually done in QFTs on curved spaces [41]. 

A simpler regularization was proposed by Liu and Mezei [10] for a spherical or any 
scalable entangling region on a flat space-time. They define the renormalized entanglement 
entropy (REE) by acting a differential operator of the radius of the sphere on the EE. In three 
dimensions, the REE of a disk of radius Rdisk becomes (1.6) which subtracts the UV divergence 
of the EE. Moreover it has been shown that the REE defined in this way is monotonically 
decreasing along any RG flow in three dimensions [11], known as the E-theorem [24, 25]. 

In our case we can define two types of REEs on the cylinder. First we note that the 
finite part of the EE of the cap-like region A on § 2 equals to that of a disk on M 2 if the 
theory is conformal. The finite part of the EE is minus the finite part of the § 3 free energy, 
F = — log 3 , as was shown by [23]. Thus the entropy Sa(9 o) for CFT 3 on the cylinder takes 
the form: 

„ m 27rI?sin#o _ . 

£U(#o)|cft = «--- F, (3.4) 

with a non-universal coefficient a. The first term, fixed by requiring the area law, is propor¬ 
tional to the circumference R sin do of the entangling surface. 

A straightforward generalization of (1.6) is to define the REE on the cylinder as 

F lm (R, e 0 ) = ( R8 r - l)S A (R)\e 0 , (3-5) 

where the derivative with respect to R is taken at fixed angle 0q. It is finite and becomes 
Elm = F at any RG fixed point thanks to the relation (3.4). We will see in section 4 that 
for a massive scalar field of mass m, the REE Elm (mR) is monotonically decreasing in the 
small and large mR regions as mR increases at fixed do. At mR = 0 it takes the value E sca i ar 
for a scalar field in three dimensions and decreases to 0 as mR —> 00 . We, however, were 
not able to determine the shape of Elm in the intermediate regime 1 -C mR -C 00 even 
numerically because of the finite lattice size effect. Thus we do not know if the ELM-hmction 
monotonically decreases along the entire RG flow of a massive scalar field theory. It is to be 
noted that this gives rise to a family of E-functions parametrised by do and they all interpolate 
between the UV and the IR fixed points of a massive scalar field in three dimensions, but 
none of them are stationary 4 as function of (mR) 2 at the UV fixed point as will be shown in 
the next section. 

4 The REE for a relevant perturbation of CFT is called stationary if the first derivative with respect to the 
coupling constant vanishes at a conformal fixed point. The REE of a disk (1.6) is known to be non-stationary 
[27, 28] for a massive free scalar theory. 
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The second way to renormalize the UV divergence of the entropy is to define the REE 
on the cylinder as 


T c (R,e 0 ) = {tan9 0 d do - 1)Sa(0 O )|* , (3.6) 

where the derivative with respect to 9 o is taken at fixed R. As the definition implies, J~c 
is always finite for the differential operator kills the area law divergence. Also it coincides 
with the finite part of the § 3 free energy F at a conformal fixed point. We will see in the 
next section for a free massive scalar field of mass m that J~c decreases monotonically as a 
function of ( mR) 2 at fixed 9q and it is also stationary as a function of ( mR) 2 at the UV fixed 
point. Then the REE Fq, obtained from EE on the cylinder, serves as an E-function in three 
dimensions. It decreases monotonically from the UV to the IR and is stationary at the UV 
fixed point for a massive scalar field. This is analogous to the Zamolodchikov’s c-function in 
two dimensions, at least for a massive scalar field. 

Before closing this section, we comment on the flat space limit of the cylinder EE. The 
cylinder metric (3.1) reduces to the flat space in the R — > oo and 9 —> 0 limits with r = R9 
held fixed, and the cap-like entangling region A turns into the disk of radius Rdisk = R9$. It 
follows from their definitions that the two REEs (3.5) and (3.6) lead to the REE of a disk 
(1.6) in this limit. 


4 Free massive scalar field 


We will calculate the EE of the cap-like region A on the cylinder for a free massive scalar 
field. We assume that the scalar field is conformally coupled to the background geometry in 
the massless limit, whose action takes the form of 


I = 



+ m 2 (j) 2 , 

O 


(4.1) 


where 1Z is the Ricci scalar. Applying the conformal decompactification and regarding the 
theory as a relevant perturbation of a free massless scalar theory by the mass term, the 
entropy will be expanded in the small mass limit and the leading term of order m 2 will be 
evaluated. On the other hand, the large mass expansion will be carried out following [29] 
and the order 1/m term of the entropy will be fixed for a general entangling surface. Finally 
the results in the two limits will be confirmed by the numerical calculation that shows the 
REE, Fc(9o,mR), monotonically decreases as mR becomes large. We also comment on the 
properties of the other REE, Elm> and the obstacles we encounter in calculating it numerically 
on lattice. 


4.1 Small mass expansion 

The EE is expected to have a series expansion with respect to the scalar mass in the small 
mass region. In order to fix the leading term of the expansion we are to calculate the derivative 
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of the free energy F n on the n-fold cover M n of K x § 2 


dm 2 n 


1 

2 



d 3 Xy/gG n (x,x) , 


(4.2) 


where G n (x,x ) is the coincident point Green’s function on A4 n . 

Using the conformal transformation (3.3) and the relation between the free energies (2.10), 
it is equivalent to that on X H 2 , 


dm 2 71 


u 


d 3 xy/§ e 2a< ' x ^G n (x,x) , 


(4.3) 


with the dilaton field 

e -2a(x) = _ sin2 _ 

(cos r + cos 9q cosh u) 2 + sin 2 9q sinh 2 u 


(4.4) 


There appears the coincident point Green’s function G n (x,x) in (4.3) which is independent 
of the position x due to the homogeneity of § 2 x H 2 . This comes out of the integral and we 
are left with the integral of the conformal factor. There is a UV divergence in the coincident 
point Green’s function which is canceled in the combination F n — n F \: 


d V r ~ - 1 

(F n - n Fi) = -y [G n (x,x)\ m 2 =0 -nG i(®,x)| m 2 =0 | + 0(m 2 ) 


(4.5) 


where V n is the integral of the conformal factor on § 2 x H 2 , 5 


V n = f d 3 x\f?je = 2n ir 3 sin 9o R 3 . (4.6) 

J S^xH 2 

There remains the coincident point Green’s function which can be obtained by construct¬ 
ing the eigenfunctions of the scalar field on x H 2 (see e.g. [18]). Inspecting the results in 
[18, 27] we find 


lim d n 

n—¥ 1 


G„ 


[.T, X)\ m 2—Q 


nG i(x,x)\ m2=0 


1 

32 R ' 


(4.7) 


Finally, the replica trick (1.4) yields the leading behavior of the EE of the cap-like region in 
the small mass limit: 

Sa(9q , mR) = a 27rRsm()o _ j? ealar - s i n g 0 ( mR) 2 + O (( mR ) 4 ) , (4.8) 

6 oZ 

with F scalar = —(ln2)/8 + 3C(3)/(16vr 2 ) w 0.0638 [25]. 

It is easy to see from the above expression that the REE J 7 lm(^ ) 9q), as defined in (3.5), 
is not stationary at the UV-fixed point mR = 0 for any value of 9q. We will discuss it in 
detail in section 4.3. 

5 The detail of the integral (4.6) can be found in appendix C. 
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4.2 Large mass expansion 

Although the analytic calculation of entanglement entropy is intractable even for free field 
theories if not conformal, one can expect to find a systematic expansion of the entropy for 
theories with a large mass gap m in powers of 1/m [ 8 , 29]: 


e | a 0 , c 2n+l 

= a— + Pm tv - 7S + 2^ 

n =0 


(4.9) 


Here ys is the topological entanglement entropy [ 6 , 7] depends on only the topology of the en¬ 
tangling surface E = dA and j3 is a scheme-independent constant [42, 43], while the numerical 
constant a is scheme-dependent. For example, a free massive scalar field has (3 = —1/12 and 
7 s = 0. The coefficients d/ n+1 are to be given by local integrals of functions of the extrinsic 
curvatures and its derivatives on the entangling surface because of the short-range correlation 
of order 1/m near the surface. We note that the expansion (4.9) has no proof for its validity 
in general, but is likely to hold for any entangling surface that is the disjoint union of a set 
of smooth curves without self-intersections [29, 44]. 

Turning to the coefficients c^ n+l for a free held theory, one can systematically determine 
c 2 n+i by the coefficient of the logarithmically divergent term of the entanglement entropy 
in the (2 n + 4)-dimensional massless free held theory compactihed on T 2n+1 [29, 45] . The 
compactihcation yields an infinite tower of massive helds in (2 + 1 ) dimensions, and the 
entanglement entropy across the entangling surface H 2 n +2 = E xT 2n+1 should equal to the sum 
of the entropies for the massive helds across E. One hnds that the sum of the entropy of the 
order l/m 2n+1 terms over the Kaluza-Klein modes gives rise to a logarithmic UV divergence, 
which should be equated with the conformal anomaly term J §( 2n+4 )| = loge in the 

higher-dimensional theory. Inspections of these coefficients lead to the following relation [29] 


s _ vr(27r) n (2n - 1)!! j2n+4) 

c 2n+l — Vol(T 2n+1 ) ' Ss 2 n +2 ’ 


(4.10) 


which can be used to determine the coefficients c^ n+1 from the conformal anomaly term in 
(2 n + 4) dimensions. 

To work it out explicitly for cf (i.e., n = 0), we start with a four-dimensional theory 
wrapped on S 1 and consider an entangling surface S 2 = S x S 1 wrapped on the S 1 . The 
entanglement entropy has a logarithmic divergence Syfili ou = s^^loge whose coefficient is 
known as Solodukhin’s formula [46, 47]: 

4 3 a +1) = |xp2! + /)/, (K^-K abab -j + k;x , '-\(K'‘) 2 ) . ( 4 -H) 

where x[^ 2 ] the Euler characteristic, 1Z the Ricci scalar, lZ aa = ^ v n c / 1 n < / 1 lZ a b a b = VJ lvpa n/j// J n° p i^ > a ^ 
and k/ JV = 7 ^ 7 is the extrinsic curvature for the normal vectors ra“ (o = 1 , 2 ) on E 2 
with the induced metric 7 ^ = g^ v — n“n“. The central charges a and c are normalized so 
that a real scalar held has (a, c) = (j|q, 77 ]) ■ In the present case where E is topologically a 
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circle, the entangling surface E 2 is topologically a torus with xfiy = 0, and only the second 
term (4.11) remains. For a circular E parametrized by 9 = 0(0), the timelike and spacelike 
unit normal vectors to E 2 are 


< = ( 1 , 0 , 0 , 0 ), 


rii = 


R sin 6 


\J sin 2 6 + (0'(0)) 


=( 0 , 1 , -&(</>), 0 ) 


(4.12) 


where the fourth components are in the S 1 direction on which E 2 is wrapped. The extrinsic 
curvature for the timelike normal vector n 1 vanishes due to the time translation invariance. 
A short calculation shows that 7 Z = 2 /R 2 , 1Z a a = 1 / R 2 , R-abab = 0, and 


1 


K 


- 2 = k 2 h llU - -(k 2 /M 2 - 


(2 cos 0 0 /2 + sin 0(sin 0 cos 0 — 0")) “ 


2 V> 7 2R 2 (sin 2 0 + 0' 2 ) 3 

Combining with (4.10) and (4.11), we find 




f s L 


1 


3 R 2 




(4.13) 


(4.14) 


A few comments are in order: 


• Our result (4.14) for the coefficient cf reproduces that of [29] in the R —> 00 and 0 —> 0 
limit with i?sin0 kept fixed, under which the entangling surface becomes a curve on 
M 2 . 


• We assumed that E is a single curve on a sphere so far, but this result (4.14) holds for 
any entangling surface which is a disjoint union of curves because the uplifted entangling 
surface E 2 in (3 + 1) dimensions is a disjoint union of tori whose Euler characteristics 
vanish and the integral (4.14) over E is just the sum of the integrals over all disjoint 
curves. 


In particular, a cap-like entangling region with opening angle 6 0 is defined by 0(0) = do, 
and the coefficient (4.14) takes the simple form: 


cap C 7T . 

A = --^-smfc'o 


1 cot 2 do 

- H-- 

3 2 


(4.15) 


Combined with the expansion (4.9) we find 


c* in 27tR sin 0 o tt vr (l sin 2 d 0 

S A [po,mR) = a --mi?sin0 O - txx— p . a x- 

e 6 lzOmxtsm^o \* o 


+ 0((mR)~ 3 ) , 
(4.16) 


in the large mR limit. 
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4.3 Numerical results 


We numerically calculated the EE of the cap-like region A on the cylinder lx § 2 by putting a 
free massive scalar field on the lattice. We closely follow the method of [31, 45] whose details 
are found in appendix B. 

Firstly, we check if the expansion (4.8) derived by using the conformal decompactification 
is valid in the small mass limit. We calculate the derivative of S A (mR, #o = ?r/2) with respect 
to ( mR ) 2 to avoid the UV divergence which contaminates our numerical precision. Fig. 1 
shows that the entropy has a linear slope 

d 2 S A (mR , 6 0 = vr/ 2 ) = -0.968 + O (( mR ) 2 ) , (4.17) 

o(mn) z 

in the small mass region mR <C 1. Integrating it by (mR) 2 leads that the entropy takes the 
form of 


S A (mR, 9 0 =-rr/2) = S A (0,d 0 = it/ 2)~ 0.968 (mR) 2 + 0((mR) 4 ) . (4.18) 

Reassuringly, this is consistent with the analytic expression (4.8) for 8q = 7 r /2 with —7r 3 /32 ~ 
-0.969. 


- -S A (0 0 =7tI2) 

d(mR) 2 



Figure 1 . The (mR) 2 derivative of the bare entanglement entropy S a (6q = 7 t/2 ) of the hemisphere 
A , in the small mass region mR « 1 . In taking the (mR) 2 derivative, we calculate bare entropies 
S A (0 o = 7 r/ 2 ) increasing (mR) 2 , and fit it as a function of (mR) 2 . The lattice size is taken as N = 1001. 

Next we examine the REE Fc(fnR,9 o) defined by (3.6) to inspect the dependences of 
the entropy on mR and 9q in the broader ranges. A detailed plot in the small mass region is 
shown in Fig. 2 (a). The REE starts decreasing from 0.0638 ~ Escaiar at the UV fixed point 
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mR = 0 for any 6q. Furthermore, it is stationary in the sense that the first derivative with 
respect to (mR) 2 vanishes at the UV fixed point: 

T c (mR, d 0 ) = F scalar + O (( mR ) 4 ) , (4.19) 

as predicted by the small mass expansion (4.8). This is in contrast to the REE T of a disk 
[10] which is not stationary at the UV fixed point of a free massive scalar theory [27, 28] 
though our Tc is supposed to reduce to T in the flat space limit. This difference may stem 
from the existence of the IR divergence on the flat space, which is regularized by the size of 
the sphere in the present setup. 

In the other extreme limit of the large mR region, we find Tc decays to zero monotonically 
as shown in Fig. 2 (b). Comparing with the expansion (4.16) which yields the large mass 
behavior of the REE 

**"*.»>- ■ < 420 > 

our numerical data are well-fitted by the curves given by (4.20) as mR becomes large. 


!Fc(0o, mR) 


0.0635 



1 1 

• 



( 


: • 

0.0630 



e 0 

• 

• 

0.0625 


• 

ttI4 (45°) 




• 

7tI3 (60°) 


0.0620 


• 

/r/2 (90°) 





J 


0.0615 





0.02 0.04 0.06 


!Tc{0o, mR) 



(mR) 2 



mR 


(a) 


( 6 ) 


Figure 2. The mR dependence of the REE Tc(mR,9o) in the small mass region mR <C 1 (a) 
and the large mass region mR 1 (b), with different cap angles 6 0 = 7r/4 (blue dots), n/3 (yellow 
dots) and n/2 (green dots). The lattice size is taken as N = 501. (a) The Tc starts from a value 
Tc(mR = 0) ~ 0.06385 (gray dotted line) at mR = 0 with vanishing slope with respect to (mR) 2 . 
This result reproduces the expected small mass expansion (4.19) Tc = E sca i ar + O ((mR) A ) in the 
small mass region mR <C 1, which means that Tc starts from the UV CFT value F sca i ar ~ 0.06381 at 
mR = 0 without any first order term of (mR) 2 . (b) It asymptotes to the leading term 7r/(120mi? sin 9 0 ) 
(solid lines) of the expected large mass expansion (4.20) Tc = n / (120mRsin 9 0 ) + 0 (l /(mR) 3 ) in the 
large mass region mR 1. 

The whole shapes of the REEs are depicted in Fig. 3. Clearly, the REEs are finite and 
monotonically decreasing to zero as m.R is increased for any 9q. Also it is a monotonic 
function of do for fixed mR, implying that increasing do from 0 to 7 t /2 can be regarded as an 
RG flow. The behavior of Tc(mR., do) is reminiscent of the REE of a disk on the flat space 
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[10, 27, 28] and the proof of monotonicity might proceed along the same lines as the proof of 
the E-theorem in [11], 


Tc{&o, mR) 


0.06 

0.05 

0.04 

0.03 

0.02 

0.01 


s 


2 


/ 


A 


&L 


• 

ttI4 

(45°) 

• 

tt/3 

(60°) 

• 

tt/2 

(90°) 

V 


J 


5 > '• i i i l 

4 6 8 10 {mR) 


Figure 3. The mR dependence of the REE Rc{mR 1 0 q) with different cap angles 9q = 7 t/4 (blue 
dots), 7r/3 (yellow dots) and 7 t/2 (green dots). It is monotonically decreasing for the all mR. In taking 
the 0o derivative, we calculate bare entropies Sa{9 o) increasing and fit it as a function of 6q. The 
lattice size is taken as N = 501. 

To recapitulate, all of the numerical results correctly reproduce both the small and large 
mass expansions in the previous two subsections, and shows that the REE J-q (3.6) always de¬ 
creases monotonically with both the scale mR and the cap size 6q increased. These numerical 
calculations give non-trivial checks for the conformal decompactification method introduced 
in section 2 and the Solodukhin’s formula (4.11) on a curved space because our numerical 
algorithm does not rely on any replica trick. 

5 Comments on 

In section 3, we introduced two types of REEs on the cylinder which lead to the REE of a 
disk (1.6) in the flat space limit, and we have solely dealt with the analytic and numerical 
properties of J~c defined by (3.6) so far. Here we will examine the Liu-Metzei type REE (3.5) 
using a free massive scalar field for comparison. Contrary to the agreements for J-q, as we 
will see soon below, there appears an incongruity between the large mass expansion and the 
numerical calculation. We will discuss possible resolutions to this puzzling situation in the 
end of this section. 

5.1 Analytic results in small and large mass limits 

The small mass expansion (4.8) of the EE of a free massive scalar held of mass m leads to 
the small mass behavior of the Liu-Mezei type REE 

7T 3 

^lm (mR, 9 0 ) = F sca lar - — sin0o(mR) 2 + O ((mR) 4 ) , (5.1) 
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which decreases linearly in ( mR ) 2 around the UV fixed point. It is not stationary in the 
Zamolodchikov’s sense as the REE on the flat space is not [10, 27, 28] while (5.1) is not a 
sensible expansion on the flat space because the 0(m 2 ) term diverges in the limit R —> oo 
and Oq —>• 0 with Rdo fixed. It clearly shows that the breakdown of the perturbation theory 
emanates from the IR divergence, the volume of the flat space [28, 48]. The REE on the 
cylinder, on the other hand, regularizes both the UV and IR divergences and is suited to the 
perturbative expansions. 

Similarly using (4.16), the leading term in the large mass expansion is given by 

- 120 J Rsineo (l - + 'O ((mnn . (5.2) 

that equals to (4.20) up to the do dependent constant. Thus it monotonically decreases as 
increasing mR for fixed do as well as increasing do for fixed mR. Both (4.20) and (5.2) 
precisely reduce to the flat space result in [27] as expected. Our analytic results (5.1) and 
(5.2) show the monotonic decrease of the Liu-Mezei type REE JVm in the small and large 
mass limits. 

5.2 Numerical results 

In the numerical calculation of J- lm, one can no longer use the same algorithm as for J-q 
due to two obstacles. One is that the definition (3.5) includes the derivative d R that re¬ 
quires the variation of the sphere radius R as opposed to the previous case. The other is 
that the discretization of the angle d by 5d = n/N causes the linear growth of the lat¬ 
tice spacing e = R5d = nR/N in R, and one cannot remove by the differential opera¬ 
tor (Ri.3 r — 1) the e dependence of the entanglement entropy because the area law term 
a{2'KRsm.do)/e = a(2!Vsin0o) becomes independent of R. To circumvent these obstacles, 
we employ a different regularization method; we calculate entropies by increasing both mR 
and N = irR/e simultaneously, that is, fixing their ratio mR/N = me/ir , and apply the 
differential operator mR. d( mR ^ — 1 = Ndjy — 1 on the fitted results. This prescription removes 
the dominant order O(N) area law term successfully. We checked that do, N) be¬ 

comes independent of N, namely, the REE J-"lm takes the same value for a different ratio 
mR/N = me/ir as long as mR is the same. In this way, the e dependence of the entropy is 
removed in the numerical calculation. 

The resultant J -lm correctly reproduces the expected small mass expansion (5.1) in 
the small mass region mR < 1, as shown in Fig. 4. The REE JVm starts from the value 
JxmI (mR) 2 =o = ^scalar — 0.0638 of the UV CFT (a free massless scalar field) like the REE 
Vc> but in a non-stationary way d^ mR yFiM\{mR) 2 =o / 0. 

We, however, unexpectedly found a discrepancy between the numerical result and the 
large mass expansion (5.2). The whole mR dependence of J-lm is drawn in Fig. 5. In the 
large mass region mR 1, The plot shows R\m asymptotes to the trivial IR CFT value 
JVm = 0, but it does not monotonically decrease in the large mass region. In fact, the 
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rLM(0O. (^) 2 ) 



(mR) 2 


Figure 4. The (mR) 2 dependence of the REE (mR, 6 q) with different cap angles 9 0 = 7r/6 
(light blue curve), 7t/4 (blue curve) and tt/2 (black curve). They correctly reproduce the small mass 
expansion (dotted magenta line) in the small mass region mR C 1. In taking the R derivative, we 
calculate bare entropies SA(niR, N) increasing both mR and N proportionally, and fit it as a function 
of R. 


numerically obtained REE J-lm obeys a different large mass expansion 

FiMimR, 0 Q ) = —— 71 . (l - 2sin 2 0 o ) + 0 ((mfi)' 3 ) , (5.3) 

120 mR sm Vq ' v 7 

instead of the expected large mass expansion (5.2), where the coefficient of sin 2 9q is 6 times 
as large as that of (5.2). For this equation to be true, the entanglement entropy should take 
the form of 


($0 j mR) | numerical 


2irR sin 6q 

= a - 

e 


—mR sin 8q 
6 


7r 

120 mR sin 


L ■ 2 

-Sill 

2 



+ O (( mR ) 3 ) 

(5.4) 


in the large mR limit. This is bigger than our expectation (4.16) by the amount of 


5 

6 


x 


c 7T sin 0o 
mR 


(5.5) 


This discrepancy would not affect Tq, because it would be removed by the operator (tan 9odg 0 — 
1) in the definition of T c- It would vanish in the flat space limit R —> oo and 0o —> 0 with R6q 
fixed as commented below (4.14), so it would not affect the result for Tq i n this limit. Both 
expressions (4.16) and (5.4) of Sa are consistent with the symmetry Sa(k/2 — 0q) = Sa(9q). 


5.3 Discussion 

We have encountered an interesting puzzle that the numerical calculation is not in agreement 
with the large mass expansion (5.2) of JTm, whose derivation is based on the large gap 
expansion (4.9), the dimensional reduction of the (3 + l)-dimensional free massless scalar 
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r LM ^o, mR) 



Figure 5. The mR dependence of the REE J-lm (mR, 9 0 ) with different cap angles 9o = 7r/6 (light blue 
curve), 7 t/ 4 (blue curve) and ir/2 (black curve). They obey a large mass expansion J r LM(wi?, 9q) = 
7r(l — 2 sin 2 0 o )/(12O mR sin 0 O ) + 0((mR)~ 3 ) (orange dotted curve) different from the expected one 
(5.2). 

field and the Solodukhin’s formula (4.11). On the other hand, our numerical calculation 
appears to be working correctly because the e dependence in J-"lm is completely removed and 
it exactly reproduces the small mass expansion (5.1) with the correct do dependence as shown 
in Fig. 4. 

The numerical result, if correct, implies that there is a missing contribution to the en¬ 
tropy in our derivation in section 4.2. A conceivable source of the discrepancy would be the 
postulated large gap expansion (4.9), where one could have included logarithmic terms of 
the forms log (m£z) and ^log(m^n). One can show the existence of such terms induces a 
new kind of UV divergences of the forms | log e and log 2 e in the uplifted free massless scalar 
theory in four dimensions. This yields a contradiction as such terms have not been observed 
in the free scalar theory. Thus the modification of the large gap expansion (4.9) is unlikely 
to remedy the situation. 

Although there exist the derivations for the formula based on the replica trick [47, 49] 
and holography [50], another possible resolution of the discrepancy would be that there is 
a missing term in the Solodukhin’s formula (4.11). Such a term, if exists, could depend 
on the curvature of the spacetime M or the derivative of the extrinsic curvature along the 
normal direction to the entangling surface to account for the difference between (4.16) and 
the conjectured form (5.4) in the large mass limit. Since the additional terms break conformal 
invariance we are not sure whether such terms are allowed to enter the formula. A similar 
problem was raised in a recent paper [51], where a total derivative term of the form □77, 
in the trace anomaly is argued to change the coefficient of the term proportional to k^ka” 
to the formula (4.11). The total derivative term is scheme-dependent in the sense that it 
can be removed by a counter term 77 2 to the Lagrangian [41], and is usually ignored. While 
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adding the total derivative does not solve the present puzzle, it leaves us a possibility that the 
numerical calculation uses a different regularization scheme from the one in the large mass 
expansion (5.2). In other words, the lattice regularization employed in section 5.2 might not 
respect the conformal invariance of the uplifted theory in four dimensions. 

To pin down the source of the discrepancy, it is desirable to investigate the dependence 
of the entropy on the background curvature, that remains well-understood yet. We hope to 
address this intriguing puzzle in future publication. 
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A Conformal transformation of a massive scalar field coupled to dilaton 


For a free massive scalar our prescription for the conformal decompactification can be verified 
easily. The action of the theory is given by (4.1). Let us denote the replica space free energy 
by F n . It is easier to first calculate the derivative of F n with respect to the mass parameter 
m 2 , which is given by the formula (4.2). The Green’s function satisfies the equation, 


■^2 / t\ 5 3 (x — x' 

—V 2 H-+ rrr ) G n (x, x) = v 


8 


Va 


(A.l) 


Now suppose we go to the conformally related manifold, A4 with metric g /w given by 

g,iu = e~ 2 a g^ . (A.2) 

It is easy to check that if we make the following transformation, 

&(x) ~ &( 

G n (x, x ) = e~G n (x, x )e~ 

then the Green’s function equations get transformed into 
M, but with a dilaton turned on: 


2 , 

the Green’s function equation on 


S 3 (x — x') 


(A.4) 
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This equation comes from the action 


I = \ j d 3 xy/$ 


g^du^dvd) + 


n~g) 


2 - 


+ me 




= Icft + \ f d 3 x y/]j m 2 e 2(j{x \jJ 2 . 


(A.5) 


For this action one can also compute the free energy and we can again write the derivative 
of the free energy as in (4.3). Use the transformations from (A.3), then it is easy to see that 


d 


d 


Q 2 F n = W o ■ (A.6) 

am am. 

We can work on the conformally transformed manifold as well. The free energies are 
the same modulo some additive constant independent of the mass parameter. Now let us see 
what is the interpretation of the tilde action. On the conformally transformed manifold, the 
dilaton couples to the massive scalar as 


I[h, t(x)] = Icft + ^ f d 3 x^m 2 e 2r(x) <)> 2 


(A.7) 


From the general arguments in the section 2, we saw that if we work on the conformally 
transformed manifold then a compensating background dilaton field has to be switched on 
which is given by, r(x) = a(x). With this dilaton field the above action transforms precisely 
into the action we got in (A.5). 


B Details of numerical calculations 

In this appendix, we summarize the numerical algorithm for calculating the entanglement 
entropy of the cap A on the cylinder M x § 2 for a conformally coupled free massive scalar 
field. 

B.l Angular decomposition 

The action is given by (4.1) with the Ricci scalar 1Z = 2/R 2 . We can regard this theory as a 
free massive scalar field theory 


\ [ d 3 xs/^g [g^d^bdvcb + m 2 ff 0 2 ] , 
z J RxS 2 


I = - 

2 

,2 — ,^2 , 1 


with an effective mass m 2 ff = rrr + The Hamiltonian is given as 


H= d8 

Jo Jo 


sin# („2J , ^ „(,\2 , (^ 0^) 2 , 2 

dtp ——— ( R vr + (dgcp ) H-;—9——|- m eS cp 


2 V v v ~r' sin 2 8 

where the conjugate momentum 7r = dt<\> satisfies the canonical commutation relation 


(B.l) 


(B.2) 


[4>(9, d>), n(Q', </>')] = —6(0 - 9')6(<j> - J)') , 

6(8 - 8')5((j) - (j)') 


(B.3) 


R 2 sin 8 
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The region A = {($,</>); 0 < 9 < Oq} has the rotational symmetry in the <j) direction, 
which allows us to reduce the space dimension to only the 9 direction by a following angular 
decomposition 


0 ) = / 1 ■ n ( ^7=T + E (0-n(0)sinn0 + 0„(0)cosra(/>) J 
V 7ri?sm0 y V2 “ J 

7i -(0,4>) = J. f + V (7T- n (0) sin n(j) + n n (9) cos ruj)) 

V7Tit d sm 0 y “J , 

In this angular decomposition, the Hamiltonian (B.2) 


(B.4) 


H = 


2R 


E 




7Tn(0)+ f(7Weff-R) 2 + wEt) ^n( 0 ) + 

y sm y / y v Vsm o j 


(B.5) 


and the commutation relation (B.3) becomes 

[K{0)aT n '{9')] = i5 nn ,5(6-0') 

B.2 Lattice discretization 


(B.6) 


We follow the discretization procedure [20]. The space coordinate 9 is discretized as 9j = jit/N 
(j = 1, 2,..., IV — 1) with dynamic variables 


(K, n i) 


N \J 2N 7rn ^ j ^) ’ O' — 1) > 


(B.7) 


(j/l,JV-l). 


In this discretization procedure, the Hamiltonian (B.5) 

^ oo / AT—1 AT—1 


h = w Y- £(ny 2 +£<Kyv„], 

n=—oo \ 7=1 i j= 1 

and the commutation relation (B.6) becomes 

[& n ,U*]=iS nn '6 jj ., 


(B.8) 


(B.9) 
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where iff™ is an 

l J 


zv( n ) _ &'("■) 


(iV — 1) x (N — 1) real symmetric tridiagonal matrix 
AT 2 sin 00/0 l / n 2 \ 

^ A “V^i-), O' l, JV - 1 ), 

sin tlj / 

(j = l,iV-2), 

OV 1, AT - 2), 


iff ™ 0 = 

JO 


iV 2 7T / 2 

—2 cos — + I {m eS R)- 


N 2 sin 6 * 3/2 


4 n)=0 ’ 


vr 2 y /2 sin 0 \ sin 62 ’ 
N 2 sin6*/ +1 / 2 
vr 2 yj sin 9j sin Qj+\ 

(K — j| > 1), 


(B.ll) 

(B.12) 


with a Z 2 symmetry K^-i,N-j = if?/ corresponding the parity symmetry 6 —> n — 9. This 
matrix ifl n ) is related to the correlation matrices X f ™' 1 = and = (n^Ilf,) as 

X(n) _ i(Ad n )) _1 / 2 and P = \{K' M) 1 / 2 . The size 6*0 of the subsystem A is chosen to be a 
half-integer in units of the lattice spacing, 6 $ = (r + 1/2)/N with an integer r. This choice 
corresponds to the free boundary condition in the continuum limit. In our calculation, we 
take N = O(10 2 ~ 3 ), which is sufficiently large for our purpose. 

The entanglement entropy of the disk S( 6 q) is obtained by using r x r submatrices 
xl n) = and p r {n) = (P^hKijKr as 

OO 

S(6 0 ) = S 0 + 2Y / Sn, (B.13) 

71=1 

where S n is the contribution from the n-th angular mode 


S n = tr [(C n + 1 / 2 ) log(C n + 1 / 2 ) - (C n - 1 / 2 ) log(C n - 1 / 2 )] , (B.14) 


with C n = y iff ™' 1 Pr"' 1 . In the following, we describe how to perform this infinite summation 
over n under controlled numerical errors. 


B.3 Large angular momentum 

In the large angular momentum limit n —> 00 , the correlation matrices X^ = \{K ( n )) -1 / 2 
and P- ^ = \{K l n l ) 1 / 2 approach almost diagonal matrices [29], for general symmetric tridi¬ 
agonal matrices if( n ) such as 

K jf = KjW + h(j ), 

4 n) = 0 (N-jI>2). 
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The products of the submatrices X^' 1 PlJ 1 ' 1 almost equal to 1/4 times unit matrix up to order 
1/n 8 . The nontrivial entries are at the lower-right corners 


(4“)Pi">V = j + ^-^ + 0( !/»•), 


n -1 


n u 


(Af)pW) r , r _ 1 = 0(l/n 6 ), 
(xWpW) r _ 1>r = 0(l/n 6 ), 


(B.16) 


where 


c(r) = 


t(r) 


4y / fc(r)/c(r + 1) ^y / fc(r) + y^r + 1)) 

AT 4 sin 2 0 r+ i/ 2 

7 r 4 (l/sin0 r + 1/sin0 r+ i) 2 
h(r) fe(r+l) 


2 ’ 


(B.17) 


= VW) ' V^h+I) + + k(r + 1) 


v / fe(r) + \Jk[r + 1) 2fe(r) 2fe(r + 1) ’ 

= ^ + ~ 2 " 2 cos^^ (sin(9 r + sin6» r+ i) 2 . 

The eigenvalues of the matrix C n = \JXr^ Pr"' > are l/2 + 0(l/n 8 ), except one eigenvalue 

c(r) c(r)b(r) 


1 

x + 


+ 0(l/n 8 ). 


(B.18) 


2 n 1 * n u 

Therefore, most of the eigenvalues do not contribute to the n-th entanglement entropy (B.14) 
up to order 1/n 8 and we obtain 


s„ = ^ (i - log 4T) + p log /rl + o ( i/n 8 ) 


n ’ 


n"* 


n u 


rr 


(B.19) 


This asymptotic formula is much faster than the direct calculation of (B.14). 

We perform the matrix trace calculation (B.14) for n less than some large angular mo¬ 
mentum n*, and use this asymptotic formula (B.19) for n > ?t* as long as S n (= 0(logn/n 4 )) 
is larger than the machine precision. The other higher modes are ignored. We can make the 
numerical error sufficiently small by taking n* large enough/’ 


C Integration of the conformal factor 

We can exactly perform the integral V n (4.6) of the conformal factor on § 4 x H 2 

p2nn poo p 

V n = dr du d£li (R 3 sinhtt) e - 2<T ( r : u ) , 

Jo Jo J S 1 

,27T f oo sinh u sin 2 6q 

= 2nirR 3 / dr du -—-^^— , 

Jo Jo (cos r + cos 6o cosh u) 2 + sin do sinh" u 


(C.l) 


See the last paragraph in the appendix of [44] for the detail to determine n*. 
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thanks to an integration formula 


f 2n dr in 

1 

1 

1 o (cos r + c) 2 + s 2 s 

is ) 2 ~ 1 

y/(c + is) 2 - 1 


(C.2) 


which we will apply with c = cos 0o cosh u and s = sin 0o sinh u. Because cos 0o cosh u ± 
z sin 0q sinh u = cosh(u ± z0q), the integration formula tells us that 



dr 

(cos r + cos 0o cosh u) 2 + sin 2 do sinh 2 u 


in 1 1 

sin 0o sinh u sinh(zz — i6o) sinh(« + i0o) 
27T coth u 
sinh 2 u + sin 2 0o 

(C.3) 


Plugging this result into the original integration (C.l), we finally obtain 


V„ = 4n7r 2 sin 0n R 3 


f 


ri( skill u) 


sinh 2 u + sin 2 0 q 


(C.4) 


= 2nir 3 sin 0 O R 3 
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